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Quantum thermal transport in two- quantum- dot system with Dzyaloshinskii-Moriya interaction 
(DM interaction) has been studied. The sign of thermal rectification can be controlled through 
changing the energy splitting or the DM interaction strength. The anisotropic term in the system can 
also affect the sign of rectification. Compared with other proposals [L Zhang et al, Phys. Rev. B 80, 
172301 (2009)], our model can offer larger rectification efficiency and show the potential application 
in designing the polarity-controllable thermal diode with a small size system (N=2). Moreover, 
quantum correlations of two-quantum-dots are investigated. We find that almost perfect quantum 
correlations can be obtained in the large temperature bias region, and quantum entanglement is 
more sensitive to the change of the DM interaction strength than quantum discord. 

PACS numbers: 



I. INTRODUCTION 

Thermal diode is a very important thermal devices [l]- 
03]. It has extensive application in phonon information Q. 
The heat current flow can be induced in a unidirectional 
way, and can be regarded as the basic part in controlling 
the large-scale heat flow. So far, many proposals have 
been raised to design a thermal diode by different models, 
such as Frenkel-Kontorova lattice and spin chains 0,0,0. 
In all these designs, the different lattice sizes will influ- 
ence the heat current flow significantly between the two 
reservoirs. In the spin chain model, the polarity of the 
thermal diode can be controlled by modulating the spin 
energy gap at different sites or changing the anisotropic 
interaction strength between the adjacent spins. The ef- 
ficiency of the thermal diode model can be measured by 
thermal rectification (JZ). In the six-spin-chain system, 
the thermal rectification 7Z can reach the value from -0.5 
to 0.3 by adjusting the system variables, and the different 
signs of 7Z means the change of the thermal diode model 
polarity. However, as studied in Ref. Q, it is not likely 
to control the thermal diode model polarity in the small 
system with the spin number N < 4. 

Recently, the problems of the interaction between the 
quantum-dot system and the thermal environment have 
attracted a lot of attention [9l-[l2j, for the quantum-dot 
system display its advantage in its scalability and long 
coherence time[i|. Due to the quantum tunnel effect 
in the quantum-dot system, the spin-spin interaction 
and th e sp in-orbit interaction will play an important 
role [M EDI. 

Moreover, when the material possesses the 



reversal asymmetry, the coexistence of the superexchange 
interaction and the spin-orbit interaction will induce a 
new interaction so called Dzyaloshinskii-Moriya (DM) 
inter act ion (lifOjj. It has been verified 0, HH that DM 
interaction can result in the change of the spin alignment 
in the material, and many problems of the quantum-dot 
system containin g th e DM interaction have been investi- 
gated already [Tonm,[H. 
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In this paper, we study extensively the heat transport 
problem in the two-quantum-dot system containing the 
DM interaction. We calculate the heat current in the 
steady state limit and analyze the heat current contri- 
butions of the different system states. We propose a 
polarity-controllable thermal diode in this two-quantum- 
dot system. Thermal rectification of this diode model 
is studied. We find that, in this two-quantum-dot sys- 
tem, thermal rectification 7Z can reach the value from -0.3 
to 0.35, through changing the DM interaction strength. 
Compared with the six- spin- chain proposal, thermal rec- 
tification 7Z in our proposal is larger than several cases in 
the spin-chain proposal. Based on this two-quantum-dot 
system, we can realize the reversal of thermal rectifica- 
tion in small system case. Besides, the quantum correla- 
tions are calculated for the steady state. We find that, 
compared with the quantum discord, the quantum entan- 
glement is more sensitive to the DM interaction. Perfect 
quantum entanglement and quantum discord can appear 
in the large temperature bias region, not only in the two 
zero reservoirs jlQj . 
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FIG. 1: (color online): A schematic representation of the 
two- quantum- dot system containing DM interaction. The two 
quantum dots interact with the thermal reservoirs at different 
temperatures Tl and Tr, respectively. 



II. THERMAL RECTIFICATION 
A. Model 

Consider a two- quantum- dot system and two thermal 
reservoirs as shown in Fig. HJ The left quantum dot 
interacts with the left reservoir only, and the right quan- 
tum dot interacts with the right reservoir only. The total 
Hamiltonian is [13, EH, H = H s + H B l + Hbr + Hsbl + 
Hsbr- Here, H s is the two-quantum-dot system Hamil- 
tonian, Hbv is the reservoir Hamiltonian, and Hsbv is 
the interaction between the reservoir v and the system. 
We consider the anisotropic interaction and the DM in- 
teraction between the two- quantum- dot system, and ne- 
glect the second order spin-orbit interaction. The system 
Hamiltonian isllOL 



J{l-iD)ala+ 



B + b 



B-b 



(1) 



where, the parameter x describes the anisotropy term in 
the system spin xy plane, the parameter D stands for the 
DM interaction between the two- quantum- dot system, 
and the energy gaps of the left spin and the right spin 



are 



B+b B-b 



respectively. 

In our case, the reservoir consists of a collection of 
non-interacting oscillators. The parameter uo n labels dif- 
ferent oscillator frequency. The system-reservoir inter- 
action strength g v n denotes the interaction between the 
spin v and the oscillator of frequency uo n in reservoir v. 
The reservoir Hamiltonian Hbv and the system-reservoir 
interaction Hamiltonian Hsb u can be written asfioj. 



Hbv = y^^uj n b+ u b nu 

n 

H SBv = S V B V = (<x+ + 0(E^ )? W + gP*b+ v ) 



(2) 



We want to investigate the system properties in the 
steady state limit. The eigenvalues and the correspond- 



ing eigenvectors of the system Hamiltonian are[l 
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Here, the state |0) and |1) are the spin of upper level en- 
ergy and spin of lower level energy, respectively. The 
coefficients N^, M^, £, and r\ satisfy the relation 

^ = (1 + T^SM" 172 ' and M± = (! + (^) 2 )" 1/2 . 
i = (b 2 + J 2 + (JD) 2 ) 1 ' 2 , rj = (B 2 + {Jx) 2 ) 1 ' 2 - 
We now rewrite the system components in the system- 
reservoir interaction Hamiltonian Hsbv in the new basis 
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The parameters n | is defined as, {S^l 

rij-Bb , J 2 X ICL 12 _ loL |2 _ 

I cR 12 _ v€+Bb , J 2 X \oR 12 _ I oR |2 _ y£-Bb _ f 2 X 
P4,2l - 2rj£ + 2^' P4,ll - l D 3,2l ~ 2r)£ 2^ ' 

We apply the master equation method to solve the 
system dynamics. We obtain the system time evolution 
function by the second order approximation in the inter- 
action picture [HI, l20j . 



dpn 



dt 



■i[HsBj(t),P(0)]r 
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Jo 



dr[H SBj {t),[H SBj {r),p{T)}} r 



(5) 



In our calculation, we assume the weak interaction of the 
system and the reservoirs. The Markov approximation 
is applied, and the reservoirs keep in its thermal equilib- 
rium throughout the evolution. Our goal is to see the 
system steady state behavior in the long time limit, we 
can neglect the coherence term in the system. Under 
these conditions mentioned above, we can get the Pauli 
master equation related to the system state population 
Pnn (n=l, 2, 3, 4), 



Where W nm is W nm 
the rates k>n — are, 
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and 



(7) 
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Here, n B (u>) denotes thermal distribution of the reser- 
voir v at temperature T„, i.e. n B {oj pq ) = 



BV^pqJ — e ^ P g/T IJ _ 1 - 

co n ) is the interaction 
strength between the reservoir v and the system. The 
system steady state population are, 



Pu 
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(8) 



with the basis { , ^2)5 ^3) , ^4)}- Expand the density 
matrix p s in the system basis {|00), |01), |10), |11)}, the 
system reduced density matrix p s is, 



/ n s 71+3 + n s r] ~ B 
I "33 2rj ^ "44 2rj 
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Given the system steady states obtained above, we can 
now formulate the expression of the heat current. 



B. Results 

Heat current. The symmetric expression of the heat 
current J is defined as[19j 

J = Tr[Jp] J= 1 -[H S bl,H s }+ 1 -[H s ,Hsbr] (10) 
in our case, the heat current can be written as 



(11) 



and P n = Tr B j[pn,n(t)]. 

Thermal rectification control by DM interaction. Un- 
der different temperature bias (AT = Tl — Tr) case, the 
control of the thermal rectification with system DM in- 
teraction strength is shown in the Fig. [21(a). The mean 
temperature of two reservoirs is, TM = 1/2 (Tl + Tr). 
As defined by Ref. [8] , thermal rectification 1Z satisfies the 
relation, 



K=(J+- J_)/max(J + , J_) 



(12) 



where the forward heat current J + is the heat current 
(from left to right) when the reservoir at higher temper- 
ature is connected to the left quantum dot, and the back- 
ward current J_ is the heat cuurent (from right to left) 
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FIG. 2: Thermal rectification (a) Thermal rectification with 
the DM interaction strength. The sign changes at critical 
value Deri. Dashed, solid and dot-dashed correspond to 
|A|T = 0.4, |A|T = 1 and |A|T = 1.6. The mean tem- 
perature of two reservoirs is TM — 1. (b) The population 
of four system eigenstates with the change of DM interaction 
strength, (c) Values of J+, n and J_, n with DM interaction 
strength. In (b) and (c), solid lines are for the case AT > 0; 
(i.e., the temperature of left reservoir is higher than the right 
one); dashed line are for AT < 0. The red, blue, green and 
brown lines correspond to the state |£2), \s&) and \sa) re- 
spectively, (d) Values of J ne t,n- The red dotted, blue dashed, 
green dot-dashed and brown solid line correspond to the state 
\si), \e 2 ), \e 3 ), and |e 4 ) respectively. AT = 0.4 and TM = 1 
in (b), (c) and (d). Other parameters are: % = 0.3, B=4, 
b=2, 3=1, F L = 1, and F R = 0.25. 



when the left quantum dot is in contact with the reser- 
voir at lower temperature. The energy gaps of the two 



quantum dots are not equal. From Fig. [2j (a), the ther- 
mal rectification 7Z will change from the positive value 
to the negative value after the DM interaction strength 
reaches a critical value D cr \. This means one can obtain 
rather different heat current through adjusting the DM 
interaction strength. The different system states popula- 
tion with the change of the DM interaction strength are 
shown in Fig. El (b). When the DM interaction strength 
is not strong, the system populates mainly at the low- 
est eigenstate \e±). The population of the eigenstate \e2) 
rises with the DM interaction strength. Values of J +?n 
and J_ 5U are shown in Fig. [2j (c). To illustrate the con- 
tribution of four different system states to the thermal 
rectification, Define J net as, 



Jn 



= j. 



+ ,n 



J- 



(13) 



where J+ ?n is the heat current contributed by the system 
state \e n ) (from left to right) when the left quantum dot 
is in contact with the reservoir at higher temperature, 
and J- y7l is the heat current contributed by the system 
state \e n ) (from right to left) when the left quantum dot 
is connected with the reservoir at lower temperature. We 
display the values of J net in Fig. [2j (d), in the region 
of the weak DM interaction strength, state ^2) provides 
the negative net heat current, and the net heat current 
contributions from state ^4) is positive. The heat cur- 
rent contribution from state ^4) is more notable than 
that from state l^)- The heat current contributed by 
state \e\ ) and \es) are negligence. So according to the 
formula of the thermal rectification above, the thermal 
rectification value is positive when the DM interaction 
strength is small. After the DM interaction strength 
reaches the critical value D cr i, the contribution of ^2) 
Jnet, 2 is positive, and the heat current of ^4) Jnet, 4 is 
negative. The major contribution is from ^4), and the 
sign of thermal rectification will change to the negative 
sign. In our two-quantum-dot system, the lowest system 
eigenstate changes from state ^4) to state ^2) when the 
DM interaction strength reaches the critical value D cr 2. 
We find that the net heat current contribution of state 
I £2) is negative, and state ^4) provides a positive con- 
tribution. When DM interaction exceeds critical value 
D cr 2, Jnet is mainly contributed by state l^)- From the 
analysis above, we find that the lowest eigenstate of the 
two-quantum-dot system makes the major contribution 
on the thermal rectification. Through modulating the 
DM interaction strength, we can change the sign of ther- 
mal rectification to control the polarity of thermal diode. 

Magnetic field effect. In our model, the magnetic field 
affect the energy gap of the quantum dot. We choose 
the parameters that the energy gap of the left dot is not 
zero, and the energy gap of the right one is zero. Different 
anisotropic system conditions are shown in Fig. 03 The 
two different situations are described in Fig. [3j (a) and 
Fig. El (b) with the DM interaction strength being zero 
and four respectively. From the figure, we find that the 
sign of the thermal rectification changes around certain 
value of external magnetic field. In different anisotropic 




FIG. 3: (color online): The relation between thermal rec- 
tification and the magnetic field strength in the four differ- 
ent anisotropic situations. Blue solid, red dashed, green dot- 
ted and orange dot-dashed correspond to % = 0.9, x — 0-3, 
X = 0.0 and % = —0.3 respectively. By setting b=B in Eq.(l), 
the energy gap of the right dot is zero. The system param- 
eters: the mean temperature TM — 1, the temperature bias 
AT = 1, J=l, Y L = l and Y R = 0.25. (a) The DM interaction 
strength is 0. (b) The DM interaction strength is 4. 



situations, the sign of thermal rectification will change 
at different critical value of external magnetic field. So, 
we can realize the reversal of the thermal diode model 
polarity by adjusting the external magnetic field in our 
two-quantum-dot system, also. 



III. QUANTUM CORRELATION 

Quantum entanglement. To qualify the quantum prop- 
erties in the steady state, we use the concurrence to mea- 
sure the quantum entanglement. Given the final X states, 
we have the analytic result of the concurrence C(p AB ) 
asJUj, 

J 



C(p AB ) =2max{0,-Vl + £ 2 |p? 1 - p* 22 \ 



1 / B 2 

2 1/ (^33 + P44) 2 - ^2"(^33 - PtJ 2 , 



^|P33-P44l 



^/(pii+py 2 -S(pfi-py 2 } 



(14) 
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FIG. 4: (color online): Quantum correlations with DM inter- 
action. Solid line: quantum entanglement; dashed line: quan- 
tum discord. Under the same parameter conditions, quantum 
entanglement rises with DM interaction more rapidly than 
quantum discord does. AT 0.4, TM = 1, x = 0-3, B=4, 
b=2, J=l, T L = 1, and T R = 0.25. 

Quantum discord. Correlation in the quantum aspects 
of the bipartite system can also be measured by quantum 
discord. The system quantum discord is defined as [22]. 

D(p AB ) = I(A : B) - C cor (p) (15) 

The quantum mutual information I(A : B) has the 
form as, I(A : B) = S(p A ) + S(p B ) ~ S(pab). The ex- 
pression S(p) is the Von Neumann entropy of the destiny 
matrix p, satisfy S(p) = —Tr(p\ogp). p A and p B are 
the reduced destiny matrix of the system X state den- 
sity matrix. The Classical correlation C CO r(p) satisfies 
the equation C cor (p) = S(p A ) — mine* [S(p\{Bi})], and 
the last term mmB i [S(p\{Bi})] is the condition entropy. 



Entanglement Discord 








FIG. 5: (color online): Quantum correlation with tempera- 
ture difference of two reservoirs. TR: temperature of right 
reservoir, D: DM interaction strength. Almost perfect quan- 
tum correlation can be achieved even though the temperature 
difference of two reservoirs is large. Parameters set in calcu- 
lation: T L 0.5, x 0-9, b=0, B=0.2, J=l, T L = 1 and 
T R = 0.25. 

We study the relation between the quantum correla- 
tion and the DM interaction strength. From Fig. 2J it 
is seen that quantum entanglement is more sensitive to 
the DM interaction strength than quantum discord, and 
perfect quantum correlation can be got within the small 
increment of the DM interaction strength. 



(a) State 1 (b) 




FIG. 6: (color online): Population of different system eigen- 
states |ei)(a), |e2)(b), \s3)(c) and |s4)(d) with different right 
reservoir temperatures and DM interaction strengths. The 
vertical axis labels the population. Parameters set: Tl — 0.5, 
X = 0.9, b=0, B=0.2, J=l, T L = 1 and F R = 0.25. 

Almost perfect quantum entanglement and quantum 
discord can be obtained in the large temperature bias 
region in our model, see in Fig. [5] The corresponding 
four system states population are shown in the Fig. OThe 
reason for the high value of the quantum correlations is 
that almost all system population are at the system state 

k 2 >. 



IV. CONCLUDING REMARKS 

We study the nonequilibrium two- quantum- dot system 
with the DM interaction extensively. The behavior of the 
heat current, quantum entanglement and quantum dis- 
cord have been investigated. We find that we can control 
the thermal rectification by two types of modulations, 
one is adjusting the system DM interaction strength, and 
the other is regulating the energy gap of two-quantum- 
dot system. This provides a possible new way to design a 
polarity-controllable thermal diode. Compared with the 
spin-chains model studied before this two-quantum- 
dot model can offer large thermal rectification when the 
size is small (N=2). This two-quantum-dot system can 
be used to build small size phonon devices, and has a 
potential application in the phonon information. 
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